IHEP-TH-25 



Light-cone QCD predictions for elastic ec/-scattering in the 

intermediate energy region 



Jun Cao 2 and Hui-fang Wu 1 ' 2 
1. China Center of Advanced Science and Technology (World Laboratory) 
P.O. Box 8730, Beijing, 100080, China 
2. Institute of High Energy Physics, P.O. Box 918(4), Beijing, 100039, China* 

(February 9, 2008) 



Abstract 



The contributions of helicity-flip matrix elements to the deuteron form factors 
are discussed in the light-cone frame. Normalized A(Q 2 ), B(Q 2 ), Gq(Q 2 ) and 
T20 are obtained in a simple QCD-inspired model. We find that G+_ plays 
an important role in Gq(Q 2 ). Our numerical results are consistent with the 
data in the intermediate energy region. 
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I. INTRODUCTION 



The electromagnetic properties of the deuteron have received extensive attention since 
it can be used to explore the quark and gluon degrees of freedom in the simplest nuclei. 
Many approaches are suggested to explain its electromagnetic form factors. Among them, 
the conventional meson-nucleon picture [0,0 gives a satisfying explanation in the low energy 
region, i.e. Q 2 < 1 GeV 2 . However, as mentioned by Arnold et al. ||, the experimental 
results are in sharp disagreement with the meson exchange calculations at higher momentum 
transfer (see Fig. 1). It means that quark and gluon degrees of freedom must be taken in 
account. Moreover, the success of dimensional counting rules [|J] implies that perturbative 
QCD(pQCD) may give correct predictions at large momentum transfers. PQCD predicts || 
that, in the Light-cone frame(LCF), the helicity-zero to zero matrix element Gq will be 
the dominant helicity amplitude at large Q 2 for elastic ed-scattering. Carlson and Gross 
have pointed out that LCF helicity-flip amplitudes, G+ and G\_, are suppressed by 
factors of Aqqb/Q and Aq CD /Q 2 , respectively. It is argued J7J that the dominance of Gq Q , 
thus the validity of perturbative QCD predictions, begins at Q 2 ~ 0.8 GeV 2 . Assuming the 
helicity-zero to zero dominance, A simple ratio of form factors for the deuteron is predicted 

G c :G M :G Q = (1 - ^) : 2 : -1, (1) 

with the kinematic factor rj = Q 2 /AM 2 . 

Unfortunately, only to the lowest order, over 300 000 diagrams containing six fermion 
lines connected by five gluons are required to obtain the full elastic erf-scattering amplitude 
in pure pQCD approach. Such a calculation can be found in the work of Farrar, Huleihel, 
and Zhang ||, which shows that the direct application of pQCD to deuteron form factor at 
experimentally accessible momentum transfers has still a long way to go. 

To make detailed predictions for deuteron electromagnetic form factors, we have sug- 
gested a model in the intermediate energy region. Based on the reduced nuclear amplitude 



defined by Brodsky and Hiller ||10|| , the model makes the points that a simple nuclear wave 
function can represent the data of the deuteron electromagnetic structure function A(Q 2 ) 
and shows a scenario where Gq is already dominant at Q 2 of 1 GeV 2 . Normalized Gq can be 
extracted from fitting the data of A(Q 2 ). It is pointed out [llj that the helicity-flip matrix 
element G+ can not be neglected in the expression of B(Q 2 ). In addition to that, we find 
that plays an important role in Gq(Q 2 ). Neglecting it will result in a contradiction 
with the data and the conventional meson-nucleon prediction in the low energy region. Due 
to the kinematical reason, the ratio of form factors in Eq. (1) should be modified even at 
extremely high Q 2 . We will model the helicity-flip matrix elements G+ and G~\__ in the 
intermediate energy region according to the pQCD predictions. The normalized structure 
functions A{Q 2 ) and B(Q 2 ), tensor polarization T 20 , and form factors GciQ 2 ) and Gq(Q 2 ) 
will be discussed on the basis of the normalized Gq and the above model. 

A QCD-inspired model for Gq will be reviewed in Sec. II. Under a phenomenological 
consideration, the model is applied to obtain normalized A(Q 2 ), B{Q 2 ) and T 20 in Sec. Ill, 
where Gq is also discussed. The numerical results and summary are presented in Sec. IV 
and Sec. V, respectively. 
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II. A QCD-INSPIRED MODEL 



For electron scattering on a deuteron, the Rosenbluth cross section is 



da 




Mott L 



A{Q 2 ) + B{Q 2 )t^ 



,6, 



(2) 



where A{Q 2 ) and B{Q 2 ) are determined by f[2 



MQ) = G 2 C + - V G 2 M + j?G; 
B(Q 2 ) = lv(l + v)G 2 M , 



(3) 
(4) 



with 1] = Q 2 /4M 2 , where M is the deuteron mass. The deuteron form factor Fd(Q 2 ) is 



defined as F d (Q 2 ) = sJA(Q 2 ). 

In the standard LCF, defined by |L3| q + = 0, q y = 0, and q x = Q, the above form factors 
can be obtained from the plus component of three helicity matrix elements: 



G c 
Gm 
Gn 



2p+(2r] + l) 
1 




2p+(2r/ + 1) 



(5a) 
(5b) 
(5c) 



To avoid the complicated pQCD calculation, a model for the deuteron form factor F^Q 2 ) 
has been suggested 0. The point is that the hard kernel at large Q 2 is assumed to be 
the perturbative amplitude for the six quarks to scatter from collinear to the initial two- 
nucleon configuration to collinear to the final two-nucleon configuration, where each nucleon 
has roughly equal momentum. Since the gluon is a color octet, the single-gluon exchange 
between two color-singlet nucleon is forbidden. In the lowest order, the quark-interchange 
is necessary in addition to the one-gluon exchange between two nucleons. For the binding 
energy of the deuteron is small, we can divide roughly the kernel into two parts. One 
represents the interchange of quarks and the gluon exchange between two nucleons, which 
transfer about half of the transverse momentum of the virtual photon from the struck nucleon 
to the spectator nucleon. Another part is the inner evolution of two nucleons. The first part 
leads to the reduced form factor of the deuteron and the latter leads to the form factors 
of two nucleons. A vector boson(color singlet) with an effective mass M& is introduced to 
represent the quark-interchange and the one-gluon exchange effects. 



Assuming the Gqo dominance in the structure function A(Q ), we get || 
G+(Q 2 ) = 2(2r/ + 1)F 2 (Q 2 /A) [ [dx] [dy\(j)\{x, Q)t H (x, y, Q)Mv> Q)» 



(6) 



where t# is the hard scattering amplitude and 
deuteron defined by 



is the body distribution amplitude of the 
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(7) 



1 „,2„2l 



The argument for the nucleon form factor, F^, is Q 2 /4 since, in the limit of zero binding 
energy, each nucleon must change its momentum from P/2 to (P + q)/2. Using Brodsky- 
Huang-Lepage prescription |TJJ from a harmonic oscillator wave function, A' exp[—^a 
in the rest frame, the body wave function can be written as 

*?*{y,l ± ) = Aexp 



1 1 2 | + m 2 



N 



2a 2 Ay(l-y 



(8) 



where A is determined by the normalization of the wave function and a by fitting the data 
of A(Q 2 ). A direct calculation of diagrams in Fig. 2 gives 



tii{x,y,Q) 



1 



xyQ 2 + M b 2 - (x - y) l M 2 xQ 2 



x) l M 2 ) 



(9) 



where g e R is an effective coupling constant. We have taken the nucleon mass to be half of 
the deuteron mass in Eq. (||). 



III. THE ROLE OF HELICITY-FLIP MATRIX ELEMENTS 

In the intermediate energy region, Gq Q dominates the charge form factor Gc, but not Gm 
and Gq, because the kinematic factor rj is still small. While 77 <C 5? the G+ contributions 
to both Gm and Gq are enhanced by a factor of l/y/2/rj. The G+_ contribution to Gq is 
enhanced by 1/rj. Although G+ and G+_ are suppressed for dynamical reason, they may 
contribute significantly to Gm and Gq because of the kinematic enhancement. Since Gc 
dominate A(Q 2 ) while r] is small, the Gq dominance works very well in determining A(Q 2 ). 
As for B(Q 2 ), the helicity-flip matrix element G+ must be taken into account ETI. In 
Gq(Q 2 ), both G+o and G\_ may be important. 

As is well known, pQCD predicts that G+ and G+_ are suppressed by factors Aqcb/Q 
and Aq CD /G; 2 , respectively. The QCD scale Aq CD is around 200 MeV. In order to explore 
the role of G+ and G+_, we assume the following relations phenomenologically : 

GXo = ^G+, (10a) 



G+_ = c£g+, (10b) 

where / and c are two parameters. These relations are expected to be reasonable as Q 2 ^ 
2MAqcd ~ 0.8 GeV 2 , where 2MAqcd is a scale that determines the Gq dominance 0. 
Ref. JTTJ reveals that G+ contributes significantly to Gm{Q 2 )- The G+ contribution to 
Gq(Q 2 ) was also discussed there. However, as mentioned above, G+_ should be taken into 
account in determining Gq due to the same reason as G+ - Detailed analysis shows that 
it plays an important role in the intermediate energy region. By assuming relation (|T(]) we 
have picked up the contributions of G+ and G+_ which may be comparable with that of 
Goo to Gm and Gq in the intermediate energy region. Nonleading contributions to G+ that 
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may enter at the same order of G+_ are neglected for kinematic reason. Inclusion of these 
terms will change the magnitude of / slightly. 

Substituting Eq. ([TOD into Eq. (|5c|), Gq becomes 



G 



Q 



2P+(2t] + 1) 



V 



77 2rj 



G, 



oo- 



(11) 



If G+_ is suppressed strongly, Gq(Q 2 ) will be negative while the pQCD begin to be valid. 
Since it is positive at the origin, there must be a node in the region of Q 2 < 1 GeV 2 . It is in 
sharp disagreement with the meson- nucleon picture (e.g. see Ref. 0) and the experimental 
data [nj. Another constraint on the parameters can be obtained from the data of Gm{Q 2 ), 
which reveals that Gm(Q 2 ) changes its sign at Q 2 = ~ 2GeV 2 . It turns out 

2t/ + (27?o - 1)/ - cf 2 = 0, (12) 

with r/ = QI/4M 2 . Combining Eq. flTID and Eq. (|12"D, it is shown that Gq will keep positive 
at any momentum transfers if c > 0.43. Different choices of c will produce very different 
G Q , but have little effect on A(Q 2 ) and B(Q 2 ). 



IV. NUMERICAL RESULTS 

At first, we assume the Gq dominance in A(Q 2 ), from which we get Eq. By fitting 
the data of A(Q 2 )(see Fig. 1) we obtain the parameters: M b = 0.5GeV, a = 0.21GeV and 
tteff = 9cs/^ = -l 5 - 

Then, given c, we can get the expressions of Gc, Gm, and Gq by substituting Eq. fllCf ) 
into Eqs. ([5]). As a demonstration, we will show the c = 1 case for simplicity. In this case / 
is 2?7o. 



G c 
Gm 
Gn 



G 



00 



2(2*7 + 1) 



16 8 2 



! 2 4 V 2 

1 77 

3 ' 3 77 



G 



00 



2(2t/ + 1) 
^oo 



2(277 



2(2t7 + 1)(1- 



4(2*7 + 1 ) 
r] z 



(1-^ 

V 



(13a) 
(13b) 
(13c) 



with 7/0 = Qq/AM 2 . The inclusion of helicity-flip matrix elements has little effect on A(Q 2 ) 
(see Fig. 1), but it changes the a e fj to 0.11. The normalized structure function B{Q 2 ) is 
given in Fig. 3, where we have chosen the parameter Qq = 1.85 GeV 2 (i.e. t/ = 0.13) from 
fitting the data. From Eq. (|13cT ) it is easy seen that Gq keeps positive for all momentum 
transfers for Gq is negative. For different choices of c, the normalized Gq and T20 are shown 
in Fig. 4 and Fig. 5, respectively. The corresponding magnitudes of / are shown there, too. 
At very large momentum transfers, say i] ^> rj Q , the ratio of form factors will be slightly 
modified to 



Gc '■ Gm '■ G 



<Q 



((1 + 



1/ + Icf) 



§77) : 2(1 + /) : -1. 



(14) 
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V. SUMMARY 



We have discussed the electromagnetic form factors of the deuteron in a QCD inspired 
model. Detailed kinematic analysis in LCF reveals that, provided the validity of pQCD, the 
helicity-zero to zero matrix element Gq dominates the gross structure function A(Q 2 ) in 
both the large and intermediate energy region, which is used to extract the normalized Gq 
from a simple model for A(Q 2 ). Further analysis shows that G+ and are also important 
to determine other form factors. In the present work, G^ and G+_ are modeled according 
to some pQCD predictions at high Q 2 . Normalized B(Q 2 ) is obtained, whose vanishing 
at Qq = 1.85 GeV 2 is used to determine the extent that helicity-flip matrix elements are 
suppressed to. To which extent the G+_ is suppressed will strongly effect the behavior of 
Gq in the intermediate energy region. We find that, If G\_ is suppressed strongly, Gq will 
change its sign twice and one of its nodes lies in the Q 2 < 1 GeV 2 region. It is contrary 
to the meson- nucleon picture and experimental data. If c > 0.43, Gq will keep positive for 
all momentum transfers in our model. Different choices of c have little effect on A(Q 2 ) and 
B{Q 2 ). As an example, we demonstrate a simple case by choosing c = 1, but a larger c 
shows a better fit for Gq and T 2 q. Since the momentum transfers of available data are not 
high enough to determine the value of c reliably, we just show a qualitative result. At very 
high Q 2 , the ratio of form factors, Gq '■ Gu '■ Gq, will be slightly modified. The ratio of Gc 
and Gq is the same as predicted by pQCD as Q 2 — > oo ||. It is apparent that Gm is bigger 
than that predicted by Brodsky and Hiller [7]. The kinematic factor yj*%q— l/y/2rj leads to 
the contributions of G+ to Gm in both the large and intermediate energy regions. 

By assuming relation fllPf) we have picked up the leading contributions of G+ and G~[_ 
which may be comparable with that of Gq to Gm and Gq in the intermediate energy 
region. There are corrections from the next to leading order contributions to G\ which 
have the form similar to G+_. For Gq these corrections can be neglected because of the 
lack of kinematic enhancement as G+_ has. For Gm they enter at the same order as G\_. 
Inclusion of these terms will diminish the magnitude of parameter / in the same way as 
increasing parameter c. Since G~\__ is not important in Gm itself (i.e. / is insensitive to c), 
our conclusions will not suffer from neglecting them. Precise / can not be determined with 
these terms unknown. It is shown in Fig. 4 and Fig. 5 that / is around 0.2. 

Since we obtain the normalized Gqq by assuming its dominance in A(Q 2 ) and relation fllDp 
by the validity of pQCD, our conclusions, except the ratio, are valid only in the intermediate 
energy region. It is found that Gc has a node |j| at Q 2 = 0.75 GeV 2 , which is predicted by 
conventional meson-nucleon method, too. Although our model can not predict the accurate 
position, it still shows that the node lies at somewhere Q 2 < 1 GeV 2 . At very large Q 2 , 
the contributions of hidden-color states should be taken into account to get a normalized 
Gq . To some extent, our model shows a smooth connection with pQCD predictions in the 
high energy region and with traditional nuclear physics conclusions in the low energy region. 
Thus it can be expected to unify the predictions for the deuteron form factors from the low 
energy to large energy region. The experiments at CEBAF will be crucial to build a more 
realistic model in the intermediate energy region. 
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FIGURES 



FIG. 1. Structure function A(Q 2 ) of the elastic ed-scattering in our model(the solid line), with 
Mb = 0.5 GeV and a = 0.21 GeV, The dashed line corresponds to the Paris potential calculation. 
Experimental data are taken from Ref. [^,[l~6|. The dotted line is our result with the corrections 
from G^Q and which is almost overlapping with the solid line. 

FIG. 2. The hard scattering diagrams. 

FIG. 3. Structure function B(Q 2 ) (the solid line). The dashed line corresponds to the Paris 
potential calculation. Experimental data are taken from Ref. [[l?]]. 

FIG. 4. The form factor Gq. Experimental data are taken from Ref. [15]. 

FIG. 5. The tensor polarization T20 with scattering angle 6 = 70°. The dashed line corresponds 
to the calculation with Paris potential. 



9 



1 f= 




0.5 / 1.5 2 2.5 3 3.5 4 



Q 2 (GeV 2 ) 




Fig.2 



-2 

10 w 



-3 

10 r o A Arnold et al. 




10 r 



0.5 1 1.5 2 2.5 3 3.5 

Q 2 (GeV 2 ) 



